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In light of the forthcoming high precision quasielastic electron scattering data from Jefferson Lab, it is 
timely for the various approaches to nuclear structure to make robust predictions for the associated response 
functions. With this in mind, we focus here on the longitudinal response function and the corresponding 
Coulomb sum rule for isospin-symmetric nuclear matter at various baryon densities. Using a quantum 
field-theoretic quark-level approach which preserves the symmetries of quantum chromodynamics, as well 
as exhibiting dynamical chiral symmetry breaking and quark confinement, we find a dramatic quenching 
of the Coulomb sum rule for momentum transfers \q\ > 0.5 GeV. The main driver of this effect lies in 
changes to the proton Dirac form factor induced by the nuclear medium. Such a dramatic quenching of the 
Coulomb sum rule was not seen in a recent quantum Monte Carlo calculation for carbon, suggesting that 
the Jefferson Lab data may well shed new light on the explicit role of QCD in nuclei. 
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Traditionally the nucleus is viewed as a collection of 
nucleons that interact via phenomenological potentials. 
This picture has proven successful since the establishment 
of the nuclear shell model and the interim has seen 
steady refinement, culminating today in sophisticated 
non-relativistic quantum-many-body approaches [1-5], 
A key assumption of such approaches is that the internal 
structural properties of the nucleons which comprise a 
nucleus are the same as those of free nucleons. However, 
with the realization that quantum chromodynamics (QCD) 
is the fundamental theory of the strong interaction, it is 
natural to expect that these nucleon properties are modified 
by the nuclear medium [' -9]. Understanding the validity 
of these two viewpoints remains a key challenge for 
contemporary nuclear physics. Should it turn out that 
nucleon properties are significantly modified by the nuclear 
medium, this would represent a new paradigm for nuclear 
physics and help build a bridge between QCD and nuclei. 
On the other hand, if the bound nucleon is unchanged this 
would shed light on colour confinement in QCD and force 
a rethink of numerous approaches to hadron structure. 

Experimental evidence for explicit quark and gluon ef¬ 
fects in nuclei remains elusive and to date the most famous 
example of such evidence - albeit not incontrovertible - is 
the EMC effect [10-12], First observed in 1982 [13], the 
EMC effect refers to a quenching of the nuclear structure 
functions relative to those of a free nucleon, and demon¬ 
strates that valence quarks in a nucleus carry a smaller 
momentum fraction than those in a free nucleon. Numer¬ 
ous explanations have been proposed, ranging from nuclear 
structure [] -17] to QCD effects [18-13], however, no con¬ 
sensus has yet been reached concerning the cause of the 
EMC effect. 

Although the EMC effect is best known, the first hints 
of QCD effects in nuclei came from quasielastic electron 
scattering on nuclear targets [24-26]. The differential cross- 


section for this process - with energy transfer uj and 3- 
momentum transfer |q| - has the form [27] 
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where q~ = — Q 2 = ur — \q\ and 9 is the electron scat¬ 
tering angle. The first experiment to separately determine 
the longitudinal (Rl) and transverse (R-t) response func¬ 
tions was performed at MIT Bates on 56 Fe [24]. The results 
suggested a significant quenching of the Coulomb sum rule 
(CSR) [28] : 
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compared to the non-relativistic expectation that Sx(|qr|) 
should approach unity for \q\ much greater than the Fermi 
momentum; on the proviso that the nucleon form factors are 
not modified by the nuclear medium [ 5] (Ge p and G En in 
Eq. (2) are the free nucleon Sachs electric form factors and 
uj + excludes the elastic peak). 1 The Bates result was soon 
verified at Saclay for various nuclei [ 6], 

Since that time the experimental situation has been in 
a state of flux. For example, a reanalysis of 12 C, 40 Ca 
and 56 Fe world data performed in Refs. [29, 30] - utilizing 
an alternative prescription for the Coulomb corrections - 
concluded that there is no evidence for quenching of the 
CSR. The analysis of the Coulomb corrections in those 
works was later challenged in Refs. [31-33] and quenching 
of the CSR over a wider range of nuclei was reported in 
Ref. [34]. This situation stands to be clarified however, 


1 Because the timelike region is not accessible in elastic scattering we define 
the CSR empirically by Eq. (2). 
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as quasielastic electron scattering has been revisited in a 
comprehensive set of measurements at Jefferson Lab [35]. 

With this in mind, in this Letter we present predictions 
for the longitudinal response function and the CSR using 
a theoretical framework formulated at the quark level. A 
rigorous consequence of this approach is that the quarks 
(and in principle gluons) - which are confined inside the 
bound nucleons - feel the presence of nearby nucleons and 
therefore nucleon properties are subtly modified by the nu¬ 
clear medium [6, 7]. Our formulation of this approach is 
based on the Nambu-Jona-Lasinio (NJL) model [9, 36, 37] - 
which is a chiral effective quark theory of QCD and is read¬ 
ily motivated by QCD’s Dyson-Schwinger equations [38], 
It provides a natural explanation for the EMC effect [23, 39- 
41 ] and has also been used to make predictions for the po¬ 
larized [23, 39] and isovector EMC effects [40-42], where 
the latter provides a significant part of the explanation of 
the anomalous NuTeV measurement of sin 2 0\v [40,43]. 
The nucleon electromagnetic current takes the form 

J^(p',p) = u N (p'){Y [F lp {q 2 ) t p + F ln (q 2 ) r„] 

+ [ F 2 p(q 2 ) t p + F 2n {q 2 ) r n ] }u N (p), ( 3 ) 

where r p ( n ) = ^ (1 ± r 3 ) and all quantities are defined at 
the same baryon density. For the calculation of the free 
nucleon form factors we follow Ref. [44] in all respects, 
including a dressed current with contributions from pions 
and vector mesons. Form factors for a nucleon in the nu¬ 
clear medium are determined in a manner analogous to 
Ref. [ 44 ], except that the quark propagator that enters the 
Poincare covariant Faddeev equation and the Feynman dia¬ 
grams that give the nucleon electromagnetic current, takes 
its in-medium form: S^ 1 (k) = f: — M — Y + ie\ where 
M is the in-medium dressed quark mass and V 1 ' the mean 
vector held, which for the form factor calculation can be 
eliminated by a shift of the integration variable. The dressed 
quark mass at a particular density is determined by the min¬ 
imum of the effective potential for nuclear matter, which 
is derived from the NJL Lagrangian using hadronization 
and path integral techniques [9]. Further discussion of this 
approach can be found in Refs. [9, 21, 40]. 

Form factor results for a free proton and a proton embed¬ 
ded in isospin symmetric nuclear matter are illustrated in 
Fig. 1 . For the free neutron form factors agreement with 
experiment is much the same [ 44 ] and medium effects are 
qualitatively similar to that of the proton. 2 At nuclear matter 
saturation density ( p B = 0.16 fm~ 3 ) we find that the proton 
magnetic moment increases by 7%, whereas the neutron 
(anomalous) magnetic moment decreases by 5%, relative to 
their free values. However, the proton anomalous magnetic 
moment remains almost constant with density; or more gen- 

2 These results are associated with the current of Eq. (3) in the mean field ap¬ 

proximation and the explicit density dependent RPA corrections are included 

via the Dyson equation of Fig. 2. 




Figure 1. (Colour online) Dirac (upper panel) and Pauli (lower 
panel) form factors for a free proton and a proton embedded 
in isospin symmetric nuclear matter (NM current) with ps = 
0.16fm~ 3 . Empirical results are from Ref. [45]. 

erally F 2p (Q 2 )/2M N shows only a slight density depen¬ 
dence because changes in F 2p are compensated for by those 
in Mat. At saturation density the proton Dirac and charge 
radii each increase by about 8%, whereas the Pauli and mag¬ 
netic radii grow by half that amount. For the neutron the 
Pauli and magnetic radii have similar increases as those for 
the proton, while the charge radius decreases in magnitude 
by 4 % and the Dirac radius more than doubles, reflect¬ 
ing the sensitivity of F\ n to small effects. For the Sachs 
super-ratio relevant to the (e, e' p) reactions on 4 He we 
find [Ge p (Q 2 )/Gm p (Q 2 )\/[Geop{Q 2 )/Gmo p {Q 2 )\ < 1 
(the subscript 0 indicates free proton form factors) with a 
magnitude consistent with experiment [ 46 , 47 ]. For the 
analogous neutron super-ratio we find a result greater than 
unity, consistent with expectation [ 48 ]. 

The longitudinal response function can be expressed as 
the imaginary part of the longitudinal polarization, which 
in isospin symmetric nuclear matter takes the form [ 27 , 49 ] 

R L {u,q) = —— ImIl L (u,q). (4) 

7T Pb 

The full result for I \ L (oj. q ) is obtained by solving the 
Dyson equation illustrated in Fig. 2 - equivalent to the 
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Figure 2. (Colour online) Dyson equation for the nucleon po¬ 
larization (shaded bubble). For the longitudinal polarization the 
photon-nucleon vertex is given by the zeroth component of the 
current operator given in Eq. (3). 


random phase approximation (RPA) - where the solid lines 
are the usual finite density nucleon propagators [27] and 
the diagram after the equality is the Hartree result. The 
nucleon-nucleon interaction, depicted in the third diagram 
of Fig. 2, is mediated by a, u and p° exchange and reads 


V NN (q)=r a (q 2 ) [1 F£(q 2 )\ 
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where [O] (i).( 2 ) = ^(i) • C( 2 ) and the subscripts label nu¬ 
cleon one and two. The quantities are —i times the 
full qq reduced f-matrices in the a, uj and p channels and 
represent the propagation of a meson in-medium, including 
its couplings to the quarks [ 44 , 49 ]. The vector and tensor 
couplings of the uj and p to the nucleon read 

F^ p) (q 2 )=F tp (q 2 )±F m (q 2 ), (6) 

where i = 1,2 and for the u-nucleon coupling we use 

Fn(F) = 9<jNN \Fi p {q~) + F ln (q 2 )] . (7) 

The normalization is rigorously given by g a NN = d ^lf [ 9 ], 
which at saturation density takes the value g a NN = 1-88. 
For the longitudinal polarization the Dyson equation gives 

U L ( U ,q) = irr-SUr-8U p L , ( 8 ) 

where the RPA correlations have the form: 
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and T p ( u ) = q 2 / \q\ 2 t p ^). All polarizations in Eqs. (8) 
and (9) are the sum of proton and neutron contributions: 
11 [ = U Lp + H Ln \ and the Hartree result reads 3 

UH = F 2 x U lx + F 1x F 2X \q\ 2 I 0 \ 
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where A = p, n. The polarizations IF^ 1 , 11/T and H P L P 

are obtained from Eq. (10) by the appropriate substitutions 
of Eq. (6). The remaining functions in Eqs. (9) and (10) are 
given in appendix A.2 of Ref. [ 7]. 

Hartree and RPA results for the longitudinal response 
function are given in Fig. 3 for |q| = 0.5 and 0.8GeV. The 
results labelled nuclear matter (NM) current are obtained by 
using the nucleon electromagnetic current operator, meson- 
nucleon form factors and nucleon propagator evaluated at 
Pb = 0.16fm -3 , whereas, for the free current results the 
electromagnetic current operator is that of a free nucleon. 
For all momentum transfers |qr|, and over all associated 
energy transfers lo, we find that the longitudinal response 
function determined with an in-medium nucleon electro¬ 
magnetic current is quenched relative to the result obtained 
using a free current. For the Hartree case this result is 
straightforward to understand from Fig. 1, and Eqs. (4), (8) 
and (10). For the relevant momentum transfers the proton 
Dirac form factor contribution to Rl(u], q) dominates, and 


3 This result agrees with Ref. [27] but corrects an error in Ref. [50], 


the observed quenching is directly associated with a softer 
F\ p in-medium. For momentum transfers \q\ < 0.7GeV 
the RPA correlations induce small corrections, shifting the 
peak in IirJ. 0 J, q) to slightly larger uj, indicating a net re¬ 
pulsive nucleon-nucleon interaction. For \q\ = 0.8GeV 
correlations play almost no role and Rl(uj, q) is completely 
dominated by the Hartree result, that is, JF£ : of Eq. (8). 
We stress that the longitudinal response function is domi¬ 
nated by the single nucleon contribution [4] and the effect 
of correlations is relatively small. Our results in Fig. 3 
show good qualitative agreement with the 208 Pb data from 
Refs. [34, 51]. In a more sophisticated calculation, includ¬ 
ing, for example, the neutron excess, Delta baryon and finite 
nuclear size corrections, the quenching of Rl{uj, q) would 
persist and is therefore a robust prediction. 

Results for the CSR of Eq. (2), using a nucleon electro¬ 
magnetic current operator evaluated at three baryon densi¬ 
ties: p B = 0, 0.1, 0.16 fm -3 , are presented in Fig. 4. For 
the free current (p B = 0) we illustrate both Hartree and 
RPA results and find that for |g| > 0.7GeV correlations 
do not materially contribute to the CSR (similar results are 
found for other baryon densities). For |q| ~ lGeV the 
CSR for the free current takes the value Sl — 0.82, which 
is considerably lower than the non-relativistic expectation 
of unity. However, if we take the non-relativistic limit of 
our result we do recover the naive expectation that the CSR 
saturates at unity for |g| much greater than the Fermi mo¬ 
mentum. Therefore, at |q| ~ lGeV we find relativistic 
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Figure 3. (Colour online) Hartree and RPA results for the lon¬ 
gitudinal response function in isospin symmetric nuclear mat¬ 
ter. Results labelled free current are obtained using the free 
nucleon electromagnetic current operator, whereas the nuclear 
matter (NM) current results use the in-medium current evaluated 
at pB = 0.16fm~ 3 . The 208 Pb data at |qj = 0.5 GeV is from 
Refs. [34, 51], 


Figure 4. (Colour online) CSR determined using a nucleon elec¬ 
tromagnetic current operator at: pB = 0, 0.1, 0.16fm~ 3 ; corre¬ 
sponding to a free nucleon current; at a density typical of 12 C [53]; 
and at nuclear matter saturation density. The data for 208 Pb is 
from Refs. [34, 51] and for 12 C from Refs. [54], both without the 
relativistic correction factor of de Forest [52] [to coincide with 
Eq. (2)]. The GFMC results are taken from Ref. [ ]. 


corrections to the CSR of roughly 20%, which in general 
are not adequately described by the relativistic correction 
factor proposed by de Forest [ 52 ], Using the nuclear mat¬ 
ter (NM) current ( p B = 0.16 fm -3 ) results in a significant 
further quenching of the CSR for \q\ > 0.5 GeV. For ex¬ 
ample, at jq| ~ 1 GeV we find that the modification of the 
nucleon form factors by the nuclear medium results in an 
additional 30% reduction in the CSR. The driver of this 
effect is the medium-induced change to the proton Dirac 
form factor illustrated in Fig. 1 . Modification of the Pauli 
form factors is less important because their contributions are 
suppressed by \q\ 2 /4M^ and -?%(„)/2ATy is largely un¬ 
changed in medium. These results demonstrate that the CSR 
is a particularly sensitive measure of even slight changes in 
the nucleon form factors, because with increasing |qr| these 
effects are cumulative. 

This dramatic quenching of the CSR has also been seen 
in other calculations [ 49 , 55 ], where the internal struc¬ 
tural properties of bound nucleons are also self-consistently 
modified by the nuclear medium. Such observations are 
consistent with many experiments on various nuclei (e.g. 
Refs. [ 24 , 26 , 34 , 51 , - 58 ]), as illustrated in Fig. 4 with 

a comparison to 208 Pb data from Refs. [ 34 , 51 ]. However, 
calculations that assume an unmodified nucleon electromag¬ 
netic current [51 - 62 ], including the state-of-the-art Green 
function Monte Carlo (GFMC) result for 12 C from Ref. [ 4 ], 
consistently find at most modest quenching of the CSR. In 
Fig. 4 the GFMC calculation is contrast with our CSR re¬ 
sult, obtained using a nucleon current evaluated at a baryon 
density typical of 12 C [ 53 ], which again finds a dramatic 
quenching; the squares are the 12 C data from Refs. [ 54 ], 
which at the largest |q| cannot distinguish between the two 
results. Both these formalisms have many compelling fea¬ 
tures. For example, the GFMC approach has had success at 


describing properties of A ^ 12 nuclei [1, 4 , 5 ], whereas 
our QCD motivated formalism provides a natural explana¬ 
tion for the EMC effect [ 23 , 39 - 41 ]. This impasse over 
the CSR stands to be resolved however, by forthcoming 
quasielastic electron scattering results from Jefferson Lab, 
at high momentum transfer and on a variety of nuclear tar¬ 
gets [ 35 ]. Verification, or otherwise, of the quenching of 
the CSR therefore promises to soon reveal critical aspects 
of the explicit role played by QCD in nuclei. 
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